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Using a dimensional reduction formula for the lattice fermion determinant we study canonical 
determinants on quenched SU(3) gauge configurations. The canonical determinants decribe a 
fixed quark number and we analyze their properties below and above the transition temperature. 
We find that above Tc the signatures of center symmetry breaking are very strongly manifest in the 
distribution of the canonical determinants in the complex plane, and we discuss possible physical 
implications of this finding. We furthermore analyze the relative weight of the different quark 
sectors below and above the transition temperature. 
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1. Introductory comments about canonical determinants 

Canonical fermion determinants det[D](^', which describe a fixed number q of quarks, are con- 
ceptionally interesting objects. On the lattice they may be obtained from the usual grand canonical 
determinant det[D(/x)] through a Fourier transform with respect to imaginary chemical potential |i : 

det[D]^'''^ = -3- r d(pe-'i'P dtt[D{pL = i(p/P)]. (1.1) 

271 J —n 

Here j3 is the inverse temperature which is given by the (periodic) temporal extent of the lattice, (p 
is the angle that parameterizes the imaginary chemical potential pL = i(p/^. The individual canon- 
ical determinants det[D]('''' appear as coefficients in the fugacity expansion of the grand canonical 
determinant det[D(/x)], 

det[D(Ai)] = ^e^^'^ det[D]('?) . (1.2) 
q 

Thus the representation with the canonical determinants is equivalent to the grand canonical for- 
mulation. In recent years several numerical simulations in the canonical formalism may be found 
in the literature [[l]] - [Q], and were reviewed at this conference |j5|]. 

Canonical determinants do not only provide an alternative approach to lattice simulations with 
finite density, but also have interesting physical properties ||6|]. The grand canonical determinant 
det[D(;U)] is a gauge invariant object and thus is a sum of products of closed loops which are 
dressed with link variables. The chemical potential is equivalent to a temporal fermionic boundary 
condition exp(/<p). This phase at the boundary is seen by the loops that wind around the compact 
time direction according to their total winding number and gives rise to a phase factor exp(/(pfc) 



for a loop that winds fc-times. The Fourier integral ( |1.1| ) over this boundary condition projects the 
grand canonical determinant to only those loops which have a net winding number ofk = q (see, 
e.g., for a more detailed discussion of these relations). The fugacity expansion (L2) thus may 



also be viewed as an expansion in terms of winding numbers of loops. 

In a similar way all gauge invariant objects may be decomposed into sectors of loops with 
fixed winding number. This has been discussed for the chiral condensate, where it is hoped that 
the "dual chiral condensate", defined as the sector of the chiral condensate with winding 1, might 
help to understand a possible relation between chiral symmetry restoration and deconfinement 



2. Center symmetry 

The canonical determinants have simple transformation properties under center transforma- 
tions, where all temporal gauge links ?74(x,fo) at a fixed time argument are multiplied with an 
element z of the center of the gauge group, i.e., Ui^{x,tQ) —>■ zU4{x,to). For the case of gauge group 
SU(3) the z are the phases z = l,exp(ib/27r/3). The gauge action and the gauge measure are invari- 
ant under center transformations. Consequently the expectation functional for the evaluation 
of observables in pure gauge theory is invariant under center transformations, as long as the cen- 
ter symmetry is not broken spontaneously. However, such a spontaneous breaking of the center 



symmetry takes place at high temperatures [11] 



2 



Determinants with fixed quark number 



Christof Gattringer 





1 1 1 1 


1 1 1 1 1 1 1 1 1 

T/T = 


1 1 1 1 

0.70 : 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

: T/T = 1 .43 : 




0.5 




C 








0.5 


0.0 










• 


0.0 
















0.5 












-0.5 




", 1 , , 


1 1 1 1 1 1 1 1 1 


, , , 1 " 




1 r 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 






-0.5 


0.0 0.5 


1.0 


-0.5 0.0 0.5 1.0 





Figure 1: Scatter plots of Polyakov loop values in the complex plane for low (Ihs.) and high temperature. 



Observables may be classified with respect to their symmetry properties under center trans- 
formations. A simple example is the Polyakov loop P which is the trace over a temporal gauge 
transporter that winds in a straight line once around compact time. As it winds once, it sees ex- 
actly one of the link variables U^{x,tQ) which are transformed with the center element z. and we 
conclude that P transforms as P — > zP- Since it transforms non-trivially, the Polyakov may be used 
as an order parameter for the breaking of center symmetry. Below the critical temperature Tc its 
expectation value vanishes, while above Tc this expectation value is finite. 

The behavior of the Polyakov loop is illustrated in Fig. |l| where we show scatter plots of the 
Polyakov loop in the complex plane. The data are from 500 quenched gauge configurations on 



8x4 lattices generated with the Liischer-Weisz gauge action 012[]. We show two ensembles with 



temperatures of T/Tc = 0.7 (Ihs. plot) and T /Tc = 1.43 (rhs.) according to the scale setting [ |T3| ] 
with the Sommer parameter. Obviously for low temperatures the values of P are compatible with 
zero, while they are non-vanishing above Tc, where the center symmetry is broken. We note that a 
true spontaneous breaking can happen only at infinite spatial volume. In that case the system will 
spontaneously select one of the three "islands" in the complex plane and only populate this one 
island. 

Similar to the Polyakov loop we also may obtain the transformation properties of our canonical 
determinants det[D]('?). We have already observed that they consist of loops with a net winding 
number of q around compact time. Thus they have a net number of q crossings of the time shce to 
where the center transformation acts and consequently transform as 



det[D 



(9) 



z'^ det[D 



det[D](^) , 



(2.1) 



where in the last step we have used that z = l,exp(ib/27r/3). 

The canonical partition sums Z^*' for a fixed quark number q are obtained as the expectation 
values of the canonical determinants using the pure gauge theory expectation functional (..)g- As 
long as the center symmetry is unbroken, we find 



Z('?) = (det[D](«)) 



c.u. 



(det[D]('?)) 



Z^"?) = for ^mod3 / 0. (2.2) 



Thus in the low temperature phase, where the center symmetry is unbroken, the canonical partition 
sums Z^^) are non- vanishing only for quark sectors with vanishing triality, i.e., when <7 is a multiple 
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Figure 2: Scatter plots of the canonical determinants det[D] ^''\q = 0,1,2 (left to right) in the complex plane. 
The top row is for T /T^ = 0.70, the bottom for T /T^ = 1 .43. 



of 3. The fact that the center symmetry is unbroken is crucial for the argument in (22), and we 
marked the step where we use that property by "c.u. !" for "center unbroken". We remark that 
the transformation properties (lA) may be combined with the center properties of observables to 
derive selection rules for observables in the center symmetric phase [^]. One finds that the total 
triality of an observable multiplied with a canonical determinant has to vanish for a non- vanishing 
contribution in the center symmetric phase. 

The situation is different in the deconfined high temperature phase where center symmetry 
is broken spontaneously. The argument used in ( ^^ ) may no longer be applied, and Z^^) can be 
non- vanishing also for q mod 3 7^ 0. In other words, above Tc also canonical determinants det[D](^) 
with non-vanishing triality can have a non-vanishing expectation value. 

In Fig. ^ we study the behavior of the canonical determinants det [D] below (top row of plots) 
and above Tc (bottom) for quark numbers q = 0,1,2 (plots from left to right). We show scatter plots 
of the values of the canonical determinants in the complex plane. The canonical determinants were 
evaluated as described in [0], using a dimensional reduction formula for determinants [14]. The 
lattice volume is 8^^ X 4, the bare quark mass parameter in the fermion determinant is set to m = 
100 MeV, and the statistics is 500 configurations for T < (top row) and 1000 configurations for 
T >Tc (bottom). 

For the zero triality case ^ = (left column of plots) the values of det[D](''^ fall on the positive 
half of the real axis and thus give rise to a positive expectation value Z^^^ both below and above 
Tc. For the non- vanishing triality sectors with q = 1,2 (center and rhs. plots) the properties below 
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and above Tc are drastically different. Below Tc the values of det[D]^'^\q = 1,2, scatter around the 
origin in a spherically symmetrical distribution. Above Tc the distribution is rather different and 
we observe the center symmetry pattern familiar from the Polyakov loop as shown in Fig. |l[ We 
stress that for the canonical determinants the pattern is even much cleaner than for the Polyakov 
loop (Both, Fig. |l| and Fig. ^ were made with the same 8^ x 4 ensembles.). 

Let us at this point address again the role of the infinite volume, which is necessary for a 
spontaneous breaking of the center symmetry. The high temperature data in Figs. |l| and ^ show 
the star-like pattern characteristic for the center broken phase. However, all three center sectors are 
populated equally and a naive averaging over all points in the scatter plots would give a vanishing 
expectation value for the Polyakov loop P, as well as the canonical determinants det[D]^'^\q =1,2 
at all temperatures - a truly uninteresting outcome. To obtain the physically relevant result one 
must take into account that in the limit of infinite spatial volume, the system selects spontaneously 
only one of the three sectors, with the other two remaining empty. In a simulation on a finite lattice 
this may be taken into account by considering above Tc the absolute values of symmetry breaking 
observables, i.e., {\P\)g and (| det[D](^) |)g. 

3. Distribution of the quark sectors 

Let us now study in more detail how the canonical determinants behave in the high temperature 
phase. In particular how sectors with different quark numbers behave relative to each other. As 
already discussed in the last section, above Tc, where the center symmetry is broken, also the 
sectors with q mod 3 7^ can have non-vanishing expectation values. As also addressed there, one 
has to average the absolute value of the determinants, as long as one is on a finite volume where the 
center symmetry cannot be broken and all three center sectors are populated equally (if a proper 
Monte Carlo update is used). 

In Fig. H we show the distribution ( | det [D] (?) | ) g / (det [D] (°) ) g as a function of q, i.e., we nor- 
malize with respect to the trivial sector with ^ = 0. The data are for lattice size 8^ x 4, a quark 
mass of m = 100 MeV, and a statistics of 100 configurations. We show the results for different 
temperatures, ranging from T = OJOTc to T = 1.43rc 

Below Tc the distribution of the absolute value of the canonical determinants shows a Gaussian 
type of behavior with a rather narrow width. This width is increasing with temperature. Above Tc 
the distribution remains Gaussian, but the width does not seem to grow any longer with temperature. 
We remark, however, at this point, that we work with a lattice size fixed in lattice units (8^ x 4), 
and change the temperature by varying the gauge coupling. Thus increasing the temperature also 
shrinks the spatial volume. This effect could mask a further widening of the distribution above Tc, 
but even if such widening persists, it is much smaller than the effect seen below Tc. A detailed 
finite volume analysis of the quark distribution must be left for future studies. 

We conclude with discussing an important consistency check: Once the canonical determi- 
nants are known, one can try to sum up the fugacity expansion and compare this sum to the grand 
canonical fermion determinant. The plots in Fig. |3| show that the canonical determinants quickly 
approach zero as \q\ increases, and a truncation of the fugacity expansion seems justified. We im- 
plemented such a test and summed the fugacity expansion, taking into account terms with \q\ up to 
typically values of 30-50. For moderate chemical potential of up to a/i ~ 0.1 we found excellent 
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Figure 3: The distribution of (| det[D]('') |)G/(det[D] W)g as a function of the quark number q for different 
values of the temperature. 

agreement (relative error less than 10^^) between the fugacity sum and the grand canonical result, 
showing that the determination of the lowest canonical determinants is sufficiently accurate. For 
larger values of the chemical potential higher terms start to contribute which would have to be 
evaluated with higher accuracy. 
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